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g ! Abstract 

j^ ^ We analyze the exact formulae for the resistance between two arbitrary notes in a rectangular 

-)— > . 

"^ ■ network of resistors under free, periodic and cylindrical boundary conditions obtained by Wu [J. 

-)— > ■ 

rH . Phys. A 37, 6653 (2004)]. Based on such expression, we then apply the algorithm of Ivashkevich, 

F^ i Izmailian and Hu [J. Phys. A 35, 5543 (2002)] to derive the exact asymptotic expansions of the 

Q ' resistance between two maximum separated nodes on an M x A^ rectangular network of resistors 

o: 

with resistors r and s in the two spatial directions. Our results is -RMxN{r,s) = c{p) ln5 + 
^ . co{p,i) + ^^1 ^2i^ with S = MN, p = r/s and C = M/N. The aU coefficients in this 

^r^ ' expansion are expressed through analytical functions. We have introduced the effective aspect 



in 

o 
o 



ratio S,eff = ^/p ^ for free and periodic boundary conditions and ^e// = y/p C/2 for cylindrical 
boundary condition and show that all finite size correction terms are invariant under transformation 



K^ , PACS numbers: 05.50.+q, 05.60. Cd, 02.30.Mv 



I. INTRODUCTION 

The calculation of the resistance between arbitrary node of infinite networks of resistors 
is a well studied subject [l|-l3|. Resistor networks have been widely studied as models for 
conductivity problems and classical transport in disordered media |#^. 

Besides being a central problem in electric circuit theory, the computation of resistances 
is also relevant to a wide range of problems ranging from random walks (see [2] and |7|, 
and discussions below), first-passage processes [8J], to lattice Green's functions [9|. Little 
attention has been paid to finite network, even though the latter are those occurring in 



real life. Recently, Wu [lO| has revisited the two-point resistance problem and deduced a 



closed-form expression for the resistance between arbitrary two nodes for finite networks with 



resistors r and s in the two spatial directions. Later, Jafarizadeh, et.al. ll| proposed an 



algorithm for the calculation of the resistance between two arbitrary nodes in an arbitrary 



distance-regular networks. However, the exact expression obtained in [10| is in the form 
of a double summation whose mathematical and physical contents are not immediately 
apparent. Quite recently Essam and Wu based on the exact expression for the resistance 
between arbitrary two nodes for finite rectangular network obtained in [lO| has derived 
the asymptotic expansion for the corner-to-corner resistance {RMxN{f,s)) on an M x N 
rectangular resistor network under free boundary conditions. For the case M = N and 
r = s = 1 they computed the finite-size corrections to the corner-to-corner resistance up to 
order N~'^: 

, , 4, ,^ _ 0.266 070 0.534 779 ^/l 
i?^.;v(l, 1) = - log iV + 0.077 318 + ^^^^ j^^ + 0[j^. 

The computation of the asymptotic expansion of the corner-to-corner resistance (in other 
word the resistance between two maximum separated nodes) of a rectangular resistor network 
has been of interest for some time, as its value provides a lower bound to the resistance of 



compact percolation clusters in the Domany-Kinzel model of a directed percolation |13 |. 

In experiments and in numerical studies of model systems, it is essential to take into 
account finite size effects in order to extract correct infinite- volume predictions from the data. 
As soon as one has a finite system one must consider the question of boundary conditions on 
the outer surfaces or "walls" of the system. The systems under various boundary conditions 
have the same per-site free energy, internal energy, specific heat, etc, in the bulk limit, 
whereas the finite size corrections are different. To understand the effects of boundary 



conditions on finite-size scaling and finite-size corrections, it is valuable to study model 
systems. Therefore, in recent decades there have many investigations on finite-size scaling, 
finite-size corrections, and boundary effects for model systems 12|, Il4|-|20|. Of particular 
importance in such studies are exact results where the analysis can be carried out without 
numerical errors. 

In this paper we will derive the exact asymptotic expansions for resistance between two 
maximum separated nodes on the rectangular network under free, periodic and cylindrical 
boundary conditions. We will show that the exact asymptotic expansion of the resistance 
between nodes of the network for all boundary conditions can be written as 

-^RmMt. s) = c{p) \nS + coip, + E ^^^^ (1) 

p=i 

where p = r/s, S = MN is the area of the lattice and ^ = M/N is the aspect ratio. The 
all coefficients in this expansion are expressed through analytical functions. We will show 
that all finite size correction terms are invariant under transformation C, -^ l/{p ^) for free 
and periodic boundary conditions and under transformation ^ — )■ 4/(p^) for cylindrical 
boundary condition, which actually means that C,eff 

^eff = C, \fp for free and periodic b.c. (2) 

^e// = i y/p/'^ for cylindrical b.c. (3) 

can be regarded as the effective aspect ratio. 

The organization of this paper is as follows: Based on the exact expression for the resis- 
tance between arbitrary two nodes for finite rectangular network under free, periodical and 



cylindrical boundary conditions obtained in 10| we express the resistance between two most 



separated nodes in terms of Ga,i3{p, with (a, /3) = (1/2, 0) and (0, 1/2) (Sec. II). We then 
extend Ivashkevich, Izmailian and Hu algorithm [1^ to derive the exact asymptotic expan- 
sions of the resistance between two maximum separated nodes on the rectangular network 
for all boundary conditions and write down the expansion coefficients up to the second order 
(Sec. III). We also discuss our results in Sec. IV. 

II. TWO-DIMENSIONAL RESISTOR NETWORKS 

An electrical network can be regarded as a graph in which the resistance Rij is associated 
to the edge between pair of connected nodes i and j. Denote the electric potential at the 



i-th vertex by Vi and the net current flowing into the network at the i-th vertex by Jj. 
When the potential difference occurs between points i and j, the current is given by the 
Ohm's law lij = (Vi — Vj)Cij, where Cij = l/Rij is the conductance of the respective link. 
By the Kirchhoff's current law total current outflow from any point in the interior is zero, 
^ ■ lij = 0, we then flnd for the voltage 

V^ = J2V,Q^/C^ (4) 

J 

where Ci = ^ . Cij and the sum is over all nodes j which are connected to i. 

The two-point resistance has a probabilistic interpretation based on classical random 
walker walking on the network. The averaging property expressed by equation (j4]) implies 
that the voltage is a harmonic function on the interior points of the graph. This makes 
the basis for the probabilistic interpretation of the voltage |2|, |2l|-|23|. The random walk 
determined by the electrical network is defined as finite state Markov chain (for more details 
see [2]) with the transition probabilities Pij that are weighted with the conductances as 
Pij = Cij/Ci- Then, when the constant voltage is applied to the graph such that V^ = 1 and 
Vb = 0, the voltage in an interior point x is determined as the hitting probability h^ that a 
walker staring at x reaches the point a before reaching b. 

Consider a rectangular M x N network of resistors with resistances r and s on edges of 
the network in the respective horizontal and vertical directions. The closed-form expression 
for the resistance -R{MxAr}(i"i, r2) between arbitrary two nodes ri = (xi, yi) and r2 = {x2, 2/2) 
for free, periodic and cylindrical boundary conditions was obtained in |lO| . 

In what follows, we will show that the resistance R{mxN}{i^,s) between two maximum 
separated nodes of the network for all above mentioned boundary conditions can be expressed 
in terms of Ga,i3{p, M, N) only. 



RM^Nir, ^) = -^ + ^ (Go,i/2(p, M, N) + Gi/2,o(p, M, N)) , (5) 



R7>cNir, ^) = V (^O'^/^^-"' ^/2' ^/2) + GiMp, M/2, N/2) , (6) 



R'^i.Nir, ^) = V (^0'V2(P, M/2, N) + G„2fl{p. M/2, N)) , (7) 



where Ga,i3{p,C,) is given by 

G^Ap,M,N) = M ReJ^fl^"^) coth MuL^^^ 
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FIG. 1: A 6 X 4 rectangular network with free boundary conditions 
for (a, /3) 7^ (0, 0). The function uj{x) is the same for all boundary conditions and given by: 



u{x) = arcsinhy^sina; 
and function f{x) is depend on boundary conditions and given by 
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A. Two-dimensional network: free boundary conditions 

Consider a rectangular M x N network of resistors with free boundary conditions and 
with resistances r and s on edges of the network in the respective horizontal and vertical 
directions. The example of a rectangular network with M = 6, A^ = 4 is shown in Fig. 
1. The resistance between the two maximum separated nodes on the network of resistors 
with free boundary conditions is the resistance between opposite corner nodes (0, 0) and 
(M — 1, A^ — 1) of the network which is given by [12| 

r(M-l) s(iV-l) 



Tjfree 
^MxN 



V,s) 



N M 

M-l N-1 



(13) 



+ 



MN 



EE 

m=l 71=1 



[cos(g^/2) cos(0„/2) - cos (M - |) 6^ cos {N - |) 0„] ' 
r~i(l — cos6'm) + s~^(l — cos0„) 



where 6m = vrm/M, 0„ = im/N. With the help of the identity 



the Eq. (fT3!) can be rewritten in the following form 

R'^':.ir,s)^-r^-^<I^ (15) 

2 y.' ^' [cos(g^/2) cos(0 J2) - COS (M - i) g^ cos (iV - ^ 0^] ' 

(m,n)^(0,0) 

Using the fact that cos (M - i) ^^ = (-l)"^cos(^™/2) and cos (A^ - |) 0„ = 
(—1)" cos(0„/2) the Eq. flT5|l can be rewritten in the following form 

<%(n.)^-^<^lfi>-ii^ (16) 
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There are two possibilities for the restriction m + n = odd to hold, namely, m-odd, n-even 
and m-even, n-odd. Splitting the sum into two parts accordingly we obtain 

M/2-1 N/2-1 

+ W^ Y. H [f{m+l/2,n) + f{m,n + l/2)] 

m=0 n=0 

M/2-1 N/2-1 M/2-1 N/2-1 

Ar v-^ v^ 9 vrn Ar v^ v^ <, irin + 1/2) 

> > cos^ -— - — — > > cos^ 



MN ^^ ^^ N MN ^^ ^^ N 

rra=0 n=0 m=0 n=0 

where 

cos^f (1 + psin^f) 

f[m,n) = ;f5^^ „ ■ (18) 

^^ ' ^ sin^Hp + psin^^ ^ ^ 

Sums of the term cos^(x) can be carried out using the identities 

N/2-1 ,^ ^ N/2-1 

v^ oTTn N 1 v^ on(n + l2) N 

y^ cos^ = + _ \ cos^ ^ — — ^^ = — (19) 

n=0 n=0 



This yields 



^MxNVy ^) 



r{M + N) s(A^-l) 



4r 



A^ M 
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(20) 



+ Xnv E E [f{m + l/2,n) + f{m,n+l/2)] 



m=0 n=0 



^Af/2-l y^Af/2-1 

It it easy to show that f{m, N — n) = f{M — m,n) = f{m, n) and thus 

N/2-1 M/2-1 ^ N~1M-1 

\ f(m. + 1 12. n^ + f(m.. n + 1 I2M = 



Now we first express double sums X]n=o ^ X]m=o ^ /("^' '^) i^ terms of X]n=o Sm=o /l'^' ''^)' 



Af/2-lM/2-l Ar_lA./_l 
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n=0 m=0 



With the help of the identities given by Eq. (TT4|) and 
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The sums E^io [/(M/2, n + 1/2) - /(O, n + 1/2)] , ^2=0 Ifi^ + V2, iV/2) + /(m + 1/2, 0)] 
can be written as 
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Plugging Eqs. (12T|) and ( 123|) back in Eq. ( |20|) we finally obtain 

Af-17V-l 

R^Ml^Nir, ^) = -^ + v^ E E [/(^ + V2, n) + /(m, n + 1/2)] 
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The sum over m in the Eq. flM|) can be carried out using the identity 24 1 



M-1 



m=0 



JJ 4 sinh'w + sin^ ^^^ = 4 | sinh (Mw + 27r/3) | 



(25) 



Note that usingmore complicated approach the identity given by Eq. (125|) has been obtained 
previously in 12| and 25 1. 

Taking the derivative over uj from the logarithm of the left and right side of the equation 
fl2S]) we obtain 
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Thus the sum over m in the Eq. ( 1241) can be carried out as 
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(27) 



(28) 



where u{x) is given by Eq. flA4p . It is easy to see that 

smh.2uj{x) = 2y^sinxY 1 + psinx^ 



(29) 



Plugging Eqs. ( 1271) and ( l28l) back in Eq. ( I24l) we obtain that Rj^'^^j^{r, s) can be written 
in the form given by Eq. ([5]). 

B. Two-dimensional network: periodical boundary conditions 



Consider a rectangular M x N resistor network with periodic boundary conditions. Using 
a closed-form expression for the resistance between arbitrary two nodes for finite network 
given by Eq. (43) of 10| we can obtain for the resistance between nodes ri = (0, 0) and 
r2 = (M/2, N/2) of the network the following expression 

1 _ nna I Mf) -I- AT rh \ 
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Using the fact that cos (Af ^^ + A^0„) = (_i)™+" the Eq. ( I30l) can be rewritten in the 
following form 
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Splitting the sum into two parts accordingly we obtain 

M/2-1N/2-1 
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The sum over m in the Eq. f l32|) can be carried out using the identities given by Eq. fl27]) 
and 023). This yields 
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Introducing function j{x) given by Eq. (ITT]) we finally arrived to the Eq. i 



C. Two-dimensional net^vork: cylindrical boundary conditions 

Consider a rectangular M ^ N resistor network embedded on a cylinder with periodic 
boundary in the direction of M and free boundaries in the direction of N. Using Eq. (46) of 
10|, the resistance between nodes ri = (0, 0) and r2 = (M/2, A^ — 1) of the network is 

Rt.N{r,s) = -'^ (34) 
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can be rewritten in the following form 
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Following along the same lines as in the case of free boundary conditions the Eq. (136|) can 
be written as 
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The sum over m in the Eq. ( 1371) can be carried out using the identities given by Eq. ( l27|l 
and (EHD- This yields 
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Introducing function f{x) given by Eq. flT2|) we finally arrived to the Eq. ([7]). 
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III. ASYMPTOTIC EXPANSION OF i?MxTv(r,s) 



In Sec. II we have shown that the resistance between two maximum separated nodes 
on an M X A^ rectangular network of resistors with various boundary conditions can be 
expressed, in terms of the function Gi/2fi{p,M,N) and Go^i/2{p,M,N) only, (see Eqs.(j5]), 
(ED and (HD). 

Based on such results, one can use the method proposed by Ivashkevich, Izmailian, and Hu 
|14| to derive the asymptotic expansion of the Ga,/5{p, M, N) in terms so-called Kronecker's 



double series 



26|, which are directly related to elliptic 6 functions (see Appendix R|) . 



After reaching this point, one can easily write down all the terms of the exact asymp- 
totic expansion for the resistance between two maximum separated nodes of the network 
(RuxNi^jS)) using Eq. (lASp . We have found that the exact asymptotic expansion of the 
RMxN^r, s) for free, periodic and cylindrical boundary conditions can be written as Eq. ([T]). 

(i) For the free boundary conditions we obtain 
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10 



where integral Jq^^x) dx is given by Eq. (lA9p . 

Thus, the coefficients C2p(p, (p=l)2,..) in the expansion Eq. ([1]) exphcitly given by 

'^^^''' ^^ = ^^-jm^^"' [^°^'^'('^^^ + K''('Vpo] (40) 

where the differential operators il2p is given by Eq. ( IA7I1 and -^'2^ (^a/pO; -^2p ' (^a/pO 
are the Kronecker's double series which can all be expressed in terms of the elliptic 9k{iy/pC) 
{k = 2,3,4) functions only (see Appendix [D|) . 

Here we list the first few coefficients in the expansion given by Eq. ([1]) 

dp) = ^ (41) 
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To simplify the notation we have use the short hand 

0fc = e,(zro), A; = 2, 3, 4, (44) 

where tq = ^^/p. 

We have also used the following relations between derivatives of the elliptic functions 

^ln^3 = j9l + T^ln^2 and T^ln^4 = 7^3' + 7^1n^2 

OTo 4 OTo OTq 4 OTo 

Note that elliptic functions 6'2, 6'3, 64^ can be expressed through the complete elliptic integral 
of the ffist kind K = K{k) and second kind E = E{k) as 
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With the help of the identities 
(9 1 



and ^ = -e'.e'; - -O^E 

OTo 



A'^ 2"'' 



dro"-- 2 

one can express all derivatives of the elliptic functions in terms of the elliptic functions 
^2, ^3, ^4 and the complete elliptic integral of the second kind E = E{k). 

For the case M = N and r = s {C, = 1, p = 1) we reproduced the result of Essam and Wu 
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with Co = 0.07731889390945876... and C2 = 0.26607044163847837.... 

(ii) For the periodic boundary conditions we obtain 
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where integral f^fix) dx is given by Eq. (jAlOp . 
The first few coefficients in the expansion are 

dp) = ^ 

^^' 27T 



VI (, 



co(p,0 = ^ 21n- + 2Cs-lne(l + p)-21nM4 
27C \ vr 

(a 
3p-l + (l + p)ro — 



n - et) 



(50) 

(51) 

(52) 



^ ((3p - ml - et) + 7rro(i + p)etel + 4ro(i + pM - ot)§^^\ne^ 



For the case M = N and r = s (^ = 1, p = 1) we obtain: 



'■^nxnI^^ ^) 



-\nN + Co + -^ + 
vr i\^ 



(53) 



with Co = 0.2078490664166084... and C2 = 0.2660704416384784... 
(iii) And for the cylindrical boundary conditions 






cyl 



MxN 



V,s) 



Vp 



71 



/»7r 

21nA^+ / (/?(x)dx + 2Cs + 41n2-21n^2(«yp^/2)04(«Vp^/2) 



2vr^p(2p)! V S 



2c\ P 



Kr'(^v^e/2) + K^r°(^v^e/2) 



(54) 
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where integral Jq^^x) dx is given by Eq. (lAlip . 

The first few coefficients in the exact asymptotic expansion are 



c{p) = ^ (55) 

co(p,0 = — 1 21n- + 2C£;-ln^(l + p)-— arctanv/p-21nM4 I (56) 
71 \ IT 7p J 

C2{p,0 = ^ {^Sp + ^ + a + P)n^y0t - el) (57) 



72 » (5 + 3p)(^| - 0t) + ^n{l + p)0t^4^ + 4ri(l + p)(^4^ - ^2)^ In ^2 



Here we have use the short hand 

ek = 9kitn), A; = 2, 3, 4, (58) 

with Ti = ^y/p/2. 

For the case M = N and r = s (^ = 1, p = 1) we obtain: 

-^^'U(^,^) = -lniV + Co + ^ + ... (59) 

S TT iV^ 

with Co = 0.36172475911729557... and C2 = -0.28448262410676656.... 

Let us now consider the behavior of the coefficients C2k{p,0 in the asymptotic expansion 
of the resistance between two maximum separated nodes on the rectangular network under 
the Jacobi transformation (see Appendix [E]). Using Eq. flE4p and Eq. flE5l) we can easily 
check that C2k{p,C,) (for all k) are invariant under transformation tq — t- 1/tq, where 

T'o = (, \fp for free and periodic b.c. (60) 

To = ^ -\/p/2 foi^ cylindrical b.c. (61) 

Using the properties of the Q -functions and of the functions K^f (see Eq. flE4p and flE5l) ) 
we can easily check from Eq. f lA5P that the lnGQ^^(p, M, A^) have the following behavior 
under the transformation Tq — ?■ I/tq: 

Gi/2,o(p, M, N) -^ Go,i/2(p, M, N) (62) 

Equations (|5]), (Q, ([7]) and (!62|) imply that the resistance R{MxN}{r,s) between two 
maximum separated nodes of the network for all above mentioned boundary conditions is 

13 



invariant under transformation tq — ^ I/tq. This actually means that ^e// given by Eqs. ([2]) 
and ([3]) can be regarded as the effective aspect ratio. 

IV. DISCUSSION 

In Fig. 2 we plot the conventional aspect-ratio (^) dependence of the finite-size correction 
term Co{^) for the resistance between two maximum separated nodes on a M x A^ resistor 
network with the free (solid lines), toroidal (dashed lines) and cylindrical (dot-dashed lines) 
boundary conditions for several values of the factor p: (a) for p = 1 and (b) for p = 4. 
We use the logarithmic scales in the horizontal axis. The finite-size correction term co{C,) 
at first decrease until ^ = ^min'- imin = 1 for p = 1 and ^min = 1/2 for p = 4. Note 
that C,min = ^1 \fp foi' arbitrary value of p. With further increase of ^ it reverses directions, 
increasing monotonically to infinity. For large enough ^ (^ ^ 1), the finite size properties 
of the resistor network with cylindrical boundary condition and those of the torus become 
the same, which means that the boundaries along the shorter direction determine the finite 
size properties of the system; for both cylindrical boundary condition and the torus, the 
boundary condition along the shorter direction is the periodic one. For small enough ^ 
(^ <^ 1), the finite size properties of the resistor network with free boundary condition and 
those of the cylinder become the same because the boundaries along the shorter directions 
for these two are the same, that is, the free boundary condition. 

In Fig. 3 we plot the effective aspect-ratio (^e//) dependence of the finite-size correction 
term cq for the resistance between two maximum separated nodes on a M x A^ resistor 
network with the free (solid lines), toroidal (dashed lines) and cylindrical (dot-dashed lines) 
boundary conditions for several values of the factor p: (a) for p = 1 and (b) for p = 4. We 
use the logarithmic scales in the horizontal axis. We can see that finite size correction terms 
Co are invariant under transformation ,^e// ~^ l/'Ce//- The finite-size correction term Co(0 at 
first decrease until ^e// = 1 for all boundary conditions and for for arbitrary value of p. With 
further increase of ^e// it reverses directions, increasing monotonically to infinity. In Fig. 4 
we plot the p dependence of the finite-size correction term co(p) for the resistance between 
two maximum separated nodes on a M x A^ resistor network with the free (solid lines), 
toroidal (dashed lines) and cylindrical (dot-dashed lines) boundary conditions for several 
values of the aspect-ratio ^: (a) for ^ = 1 and (b) for ^ = 4. The finite-size correction term 
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log f 




log f 



FIG. 2: Conventional aspect-ratio (^) dependence of finite-size correction term cq for the resistance 
between two maximum separated nodes on a M x A^ resistor network with the free (sohd hnes), 
toroidal (dashed lines) and cylindrical (dot-dashed lines) boundary conditions: (a) for p = 1 and 
(b) for p = 4. We use the natural logarithmic scales for the horizontal axis. 




log f,: 




log f.££ 



FIG. 3: Effective aspect-ratio (^e//) dependence of finite-size correction term cq for the resistance 
between two maximum separated nodes on a M x A^ resistor network with the free (solid lines), 
toroidal (dashed lines) and cylindrical (dot-dashed lines) boundary conditions: (a) for p = 1 and 
(b) for p = 4. We use the natural logarithmic scales for the horizontal axis. 

co(p) at first decrease until p = pmin- Note that value of pmin depends on the boundary 
conditions as well on the value of the aspect ratio C,- With further increase of p it reverses 
directions, increasing monotonically to infinity. 
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FIG. 4: The (p) dependence of finite-size correction term cq for the resistance between two max- 
imum separated nodes on a M x A^ resistor network with the free (sohd hnes), toroidal (dashed 
lines) and cylindrical (dot-dashed lines) boundary conditions: (a) for ^^ = 1 and (b) for ^ = 4. 

In the present paper, we study the two-point resistor problem on planar MxN rectangular 
lattices with free, periodic and cylindrical boundary conditions. Using the exact expression 
'or the resistance between arbitrary two nodes for finite rectangular network obtained in 



10| and the IIHs algorithm [lj|, we derive the exact asymptotic expansion of the corner-to- 



corner resistance on the rectangular network for all above mentioned boundary conditions. 
All corrections to scaling are analytic. 
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Appendix A: Asymptotic expansion of Gq^(/9, 7W, AA) 
Using the expansion of the cothx 

oo 

cothx = 1 + 2^ 



-2mx 



m=l 
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we can transform the Eq. ([8]) in the following form 

M-l 

G^AP,M,Af) = ReA^5^/(z^)coth[A1/(^^)+z7r/3 



n=0 



(Al) 



A^-l 



M-l oo 



n=0 n=0 m=l 

where /(a;) is given by Eqs. (fTOl) . (fTTl) and (fT2l) for free, periodic and cylindrical b.c. 
respectively. 

Using Taylor's theorem, the asymptotic expansion of the f{x) for all boundary conditions 
can be written in the following form 



X 



i + E 



^2p 2p 



(A2) 



where K2 = p — 5/3, K4 = — 3p^ — 14p + 67/15, etc for free boundary conditions K2 = 
—p + 1/3, K4 = 9p^ + 2p + 7/15, etc. for periodic boundary conditions and K2 = —p — 5/3, 
fi;4 = 9p^ + 14p + 67/15, etc. for cylindrical boundary conditions. 
Note that function f{x) can be represent as 



1 f °° 

f{x) = - exp <^ ^ 



^2p ^2p 



(2p)!- 



(A3) 



where coefficients e2p and /€2p are related to each other through relation between moments 
and cumulants (Appendix [B]) 



K2 = £2 



K4 — S4 -\- 3£n 



We will need also the Taylor expansion of the u{x) given by Eq. IQ 



u{x) = X I x + y^ 



X 



2P ^2p 



u (^py- 



X 



(A4) 



where A = ^p, A2 = —^/pil + p)/3, A4 = y/p{l + lOp + 9p^)/5, etc. for all boundary 
conditions. 

In what follows, we shall not use the special values of these coefficients assuming the 
possibility for generalizations. 
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The asymptotic expansion of Ga,ii{p,M.,N') can be derived in the similar way as it has 
done in Ref. [ij] for the second derivative of the partition function with twisted boundary 
condition (see Eq. (18) in Ref. [ij]). Following along the same line as in Ref. IJ] we can 
obtain for the asymptotic expansion of G a^jsip , M , J^) the following expression 



GaAP^M,Af) 



2S 



vr 



lnA/'+- / (^(x)da; + CB + 21n2-21n|^^,;3(iAOI 



'0 

I ( d d^ 

+ o ^27TT + A2 



■d\ 



d\^ 



In 



Oa,pm) 



p=2 



s 



2p 



p{2p)\ 



vim 



Re Kfm) 



(A5) 



where S = A4J\f, ^ = Ai/M, Ce is the Euler constant, tj^t) is the Dedekind - rj function 



7^(r) = e"'"/i2 Y[[l- e^'^'™] 



(A6) 



n=l 



^2p+2ij) is Kronecker's double series (Appendix [D]) and 9a,/5i'T) is elliptic theta function 
(Appendix E]). 

The differential operators ^^2^ that have appeared here can be expressed via coefficients 

(^2p = £2p + ^2p-Qx ^S 



fio 



VL, 



0J2 



U4 + SUo 



(A7) 



The function if{x) is defined as 



^{x) = f{x) 



X vr — X 



(A8) 



Thus, the Jq^p{x) dx depend on the boundary conditions and given by 

2 p-1 
ip{x) dx = — l + 21n ln(l + p)h arctan^/p for free b.c, (A9) 



2 p-1 
-1 + 21n ln(l + p) -\ arctan^/p for free b.c, 

TT ^/p 



(p{x)dx = 2 In ln(l + p) 



TT 



2 2 

{p{x) dx = 2 In ln(l + p) arctauiVp 

TT ^ 



for periodic b.c, (AlO) 
for cylindrical b.c. (All) 



We are interested in the asymptotic expansion of Ga,i3{p,M.,Af) with (a,/3) = (0, |) and 

;io). 
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Appendix B: Relation bet^veen moments and cumulants 



Moments Z^ and cumulants F^ which enters the expansion of exponent 



OO h. 

X 



exp<; ^ — Fk 



OO L. 

X 



fc=l 



k\ 



l + J2:^Zk 



fe=i 



k\ 



are related to each other as 

Zi - 
Z2 - 
Z3 - 

Z, - 



27| 



F2 + Fl 

F3 + 3F1F2 + F^ 

F4 + 4F1F3 + 3F2 + QFIF2 + Fi^ 



r=l 



-^fei \ / Fk, 



2^2^ [ hj ■ ■ ■ I i- 1 



h 



k\ 



2l! . . . Zrl 



where summation is over all positive numbers {ii ■ ■ ■ v} and different positive numbers 
{ki, . . . , kr} such that kiii + . . . + kri,. = k. 

Appendix C: Elliptic Theta Functions 

In this appendix we gather all the definitions and properties of the Jacobi's 6'-functions 
and Kronecker's double series needed in this paper. We adopt the following definition of the 
elliptic 6'-functions: 

&a,i3{z, t) = ^ exp <^ irir (n + ^ - a) + 2-ni [n + \ - a) [z - \ + P)\ 



r]{T) exp {nirB^ + 2'Ki{\ -a){z~\ + P)] 



X 



CXJ 

n[i- 



,27rJr(7i+a)-27ri(z+/3) 



n=0 



1 _ 27rJT(n+l-a)+27rJ(z+/3) 



(CI) 



where B2 is the Bernoulli polynomial B2 = a^ — a + \ and ri{T) is Dedekind //-function: 



t]{t) = e 



■KiTB2/2 



n[i 



— e 



2'jTiTn 



(C2) 



n=l 



where B2 = B^ = 1/6. 



19 



The elliptic 6'-functions satisfies the heat equation 

^S-.M^.-J^i^"..!^.-) (C3) 

The relation of the functions 6a,i3{z, r) with the usual ^-functions 6i{z, r) i = 1, . . . , 4 is the 
following 

Oo,oiz,T) = 9i{z,t) 

^1,0(2;, r) = e4{z,T) 
OiA^,r) = 9s{z,t) 



2'2 



In this paper we will only need these functions evaluated at 2; = and r = ztq is a pure 
imaginary aspect ratio. To simplify the notation we will use the shorthand 

9i{iTo) = 9i{0,iTo) 
Appendix D: Kronecker's Double Series 



Kronecker's double series can be defined as 26 1 



„} „—2ni(nQ+m(3) 



-l-nif ^-^ (n + rmf 

' m.neZ ^ > 

{m,n)7^{0,0) 



y it 



In this form, however, they cannot be directly applied to our analysis. In [1^ it was shown 

that Kronecker's double series can be written as 

00 

The Kronecker's double series i^|^ (r) for the cases (a,/3) = 

(0,0), (0,1/2), (1/2,0), (1/2, 1/2) can all be expressed in terms of the elliptic ^(r) 
functions only. 

Equations for K2p (t) with p = 2, 3, 4, 5 and other useful relations for elliptic 6'-functions 



and Kronecker's double series can be found in Refs. |l4J-|l7| 
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Here we write down the Kronecker's double series K2' (iro) and -K'gp ' (^'^0) that have 
appeared in our asymptotic expansions 

IT Otq r]{tTo) 12 

n Otq r]{tTo) 12 

K|'°(^ro) = ^[|^!(^ro)-^|(zro)^|(^^o)], 

K6'^(^ro) = i-Mitro) + ei{^r,)][% el{zr,) + ^|(^ro)^|(^ro)], 

K|'°(^ro) = -^[^|(zro) + ^K^To)][fi0|(zro) + 0|(^ro)^|(zro)], 



Note that when tq — ;■ 00 we have hmits O^liro) -^ 1, 9^{iTQ) — )■ 1 and ^2(^to) — ^ 
2exp (— ^) — )■ 0, and Kronecker's double series reduce to the Bernoulli polynomials: 

hm Kf{iT,) = 5°,. 

The case tq — t- can be obtained by using Jacobi's imaginary transformation of the ^(irg) - 
functions. In this case 92{iTo) -^ ^, Osiiro) -^ -^ and 6*4(2x0) -^ ;;^exp (^-4^] -^ and 
the Kronecker's function can again be reduce to the Bernoulli polynomials. 

Appendix E: Jacobi transformation 

We also need the behavior of the 6 functions, Dedekind's r^-function and the Kronecker 
functions i^gp under the Jacobi transformation 

t^t' = -1/t. (El) 



The result for the 6 functions and Dedekind's r^-function when 2; = is given in ref. 



e,{r') = {-^rf/^e,{r), (E2) 

21 



The result for the Kronecker functions K2' and i^2p ' can be obtain from the relation be- 
tween coefficients in Laurent expansion of the Weierstrass function and Kronecker functions 
(see Appendix F in [IJ]) and is given by 

Kl'^\r') = T'^Kli''\T), (E3) 

In particular, in the case of pure imaginary aspect ratio r = iro, the 6'-functions and K2p , 
K2P ' -functions transforms under (lEip as follows 

^2(V^o) = rJ/'^4(^ro), (E4) 

ri{i/To) = TQ^^riiiTo), 
K',i'\ilT,) = {ir,f^Kll'\iT,), (E5) 
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